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ABSTRACT 

We cross-correlate QSOs from the 2dF QSO Redshift Survey with groups of galax¬ 
ies. In the southern region of the 2dF we utilise galaxies from the APM Survey. In the 
northern strip, galaxies are taken from the recent Sloan Digital Sky Survey Early Data 
Release. Both galaxy samples are limited to a depth B < 20.5. We use an objective 
clustering algorithm to detect groups in these galaxy catalogues. 

A 3(7 anti-correlation is observed between 2dF QSOs and galaxy groups, confirm¬ 
ing the effect found by Boyle, Fong & Shanks in an independent dataset. This paucity 
of faint QSOs around groups cannot be readily attributed to a selection effect and is 
not due to restrictions on the placement of 2dF fibres. By observing the colours of 
QSOs on the scales of the anti-correlation, we limit the influence intervening dust in 
galaxy groups can have on background QSO flux, finding a maximum reddening on 
the scale of the anti-correlation of E{bj — r) < 0.012 at the 95 per cent level. The 
small amount of dust inferred from the QSO colours would be insufficient to account 
for the anti-correlation, supporting the suggestion by Groom & Shanks that the signal 
is likely caused by weak gravitational lensing. The possibility remains that tailored 
dust models involving grey dust, heavy patches of dust, or a combination of dust and 
lensing, could explain the anti-correlation. 

Under the assumption that the signal is caused by lensing rather than dust, we 
measure the average velocity dispersion of a Singular Isothermal Sphere that would 
cause the anti-correlation, finding a ~ 1150 km s“^. Simple simulations reject a ~ 
600 km s“^ at the 5 per cent significance level. We also suppose the foreground mass 
distribution consists of dark matter haloes with an NFW profile and measure the 
typical mass within 1.5 h~^ Mpc of the halo centre as Mi s = 1.2±0.9 x 10^® h~^ Mq. 

Regardless of whether we utilise a Singular Isothermal Sphere or NFW dark mat¬ 
ter profile, our simple lensing model favours more mass in groups of galaxies than 
accounted for in a universe with density parameter Qm = 0.3. Detailed simulations 
and galaxy group redshift information will significantly reduce the current systematic 
uncertainties in these flm estimates. Reducing the remaining statistical uncertainty in 
this result will require larger QSO and galaxy group surveys. 


Key words: surveys - quasars, quasars: general, large-scale structure of Universe, 
cosmology: observations, gravitational lensing 


1 INTRODUCTION than evidence for non-cosmological redshifts (Bukhmastova 

2001; Benitez, Sanz & Martinez-Gonzalez 2001). Observa- 
Claims of QSO-foreground galaxy associations are now more tionally, however, the situation has been complex, with a 
usually interpreted in terms of gravitational lensing rather 
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variety of effects recorded in the literature. For example 
Williams & Irwin (1998) detected a positive correlation 
between high-redshift LBQS QSOs and APM galaxies on 
degree scales, yet Martinez et al. (1999) find no strong 
quasar-galaxy angular correlation on similar scales. Ferreras, 
Benitez & Martinez-Gonzalez (1997) found a large anti¬ 
correlation between galaxies and optically selected QSOs 
near the NGP, and suggested a difficulty in selecting QSOs in 
densely populated areas. Samples of QSOs utilised in these 
and earlier papers frequently suffered from either inhomo¬ 
geneity or a dearth of data (see, e.g., Norman & Williams 
2000, for a review). Recently, Norman & Impey (2001) have 
found a significant positive correlation between a homoge¬ 
neous sample of 90 radio-bright QSOs (with median opti¬ 
cal magnitude V ~ 18.5) and galaxies lying outside of rich 
clusters. The amplitude of this angular correlation was con¬ 
sistent with a lensing explanation. The completed 2dF QSO 
Redshift Survey contains a UVX-selected homogeneous sam¬ 
ple of around 23,000 QSOs. Contemporary models suggest 
that statistical lensing should become a stark cosmological 
effect in such vast surveys (Menard & Bartelmann 2002). 

Following Shanks et al. (1983), Boyle, Fong & 
Shanks (1988, henceforth BFS88) quoted a signfficant anti¬ 
correlation between galaxies in objectively selected clusters 
and faint UVX objects. Initially, BFS88 interpreted the anti¬ 
correlation as an effect caused by a small amount of dust in 
foreground galaxy groups obscuring background QSOs. Fer¬ 
guson (1993) and Maoz (1995) generally restricted reddening 
in clusters and rich groups, results that marginally suggested 
insufficient dust to induce the observed lack of QSOs around 
galaxy groups. This prompted Groom & Shanks (1999) to 
recast the anti-correlation signal in terms of statistical lens¬ 
ing. Rodrigues-Williams & Hogan (1994) first discussed the 
BFS88 anti-correlation result as a possible effect of gravita¬ 
tional lensing in a paper conhrming that bright QSOs are 
positively correlated with galaxy clusters. Significant posi¬ 
tive correlations between galaxies and bright QSOs, which 
cannot be explained away by dust in galaxies, continue to 
be detected (Williams & Irwin 1998; Norman & Williams 
2000; Norman & Impey 2001; Gaztanaga 2002). 

Gravitational lensing can satisfactorily explain QSO- 
galaxy associations. Interim mass lenses the area behind 
it, influencing a sample of distant objects in two related 
ways. Firstly, sources are magnified. Secondly, the apparent 
sky density of sources of given intrinsic luminosity drops. 
If the number-magnitude relation of the considered sam¬ 
ple is steep, the first effect dominates and we observe the 
more numerous, fainter QSO population. If the number- 
magnitude slope is shallow, the second effect prevails and 
we see fewer QSOs, as the area that they populate dilutes. 
Ultimately, this magnihcation bias will increase the corre¬ 
lation of QSOs and galaxies if we consider a sample near 
an intrinsically steep part of the QSO number-magnitude 
relation, and it will induce a paucity of QSOs around fore¬ 
ground mass where QSO number counts flatten (Wu 1994). 
The optical QSO number counts are steep for bright QSOs 
and flatten significantly at faint magnitudes. Observation- 
ally, then, we would expect both positive and negative cor¬ 
relations between QSOs and galaxies, depending on the ap¬ 
parent luminosity of the QSO sample. 

Gravitational lensing techniques have fast become pop¬ 
ular methods of tracing the mass distributions in our Uni¬ 


verse; see Wambsganss (1996) for a review. The attraction 
is obvious - no assumptions need be made about the virial- 
isation of the observed matter, or about how the luminous 
matter traces the underlying distribution of mass. Recon¬ 
structing individual cluster masses from the shear and mag¬ 
nihcation of background galaxies is a well-established cosmo¬ 
logical tool (Squires & Kaiser 1996 and references therein) 
and measurements of the statistical inhuence of poor groups 
on the background galaxy population can probe the distribu¬ 
tion of matter in the held (Hoekstra et al. 2001). Broadhurst, 
Taylor and Peacock (1995) have suggested using magnihca¬ 
tion bias to measure the ehect dark matter in a cluster has 
in lensing the overall distribution of background galaxies. 
Studying the inhuence clusters exert on background QSOs 
is, perhaps, an even better approach; although QSO samples 
are signihcantly smaller than galaxy samples, background 
QSOs are more readily distinguished from cluster members 
than background galaxies. 

Questions remain about the anti-correlation between 
faint UVX objects and galaxies in groups detected by 
BFS88. Is the result reproduced in a different sample? Is 
it affected by dust in galaxy groups? Is it a selection effect 
or a systematic? Is statistical lensing a viable explanation 
given the large amplitude of the anti-correlation? In this 
paper, we address these questions by measuring the two- 
point correlation function between objects in the 2dF QSO 
Redshift Survey (Groom et al. 2001) and galaxy groups de¬ 
termined from the APM catalogue (Maddox et al. 1990a) 
and from the Sloan Digital Sky Survey Early Data Release 
(Stoughton et al. 2002). The 2dF QSO redshift survey con¬ 
tains a large, homogeneous, objectively determined sample 
of objects. The QSOs are spectroscopically conhrmed, mean¬ 
ing contamination is extremely low. The 2dF survey also 
measures colour, allowing strict limits to be placed on the 
ehect of dust in galaxies on background QSOs. Finally, in 
modelling our lensing signal, we consider a diherent cluster 
profile to BFS88 and adapt our analysis so that the model 
more fairly represents the data. 

This paper primarily deals with the cross-correlation 
of QSOs versus foreground galaxy groups and what it in¬ 
dicates. In the following section we discuss the 2dF QSO 
Redshift Survey and derive necessary results from it. In Sec¬ 
tion 3 we outline the galaxy samples we cross-correlate with 
QSOs and the objective method by which we derive galaxy 
groups. In Section 4 we display the cross-correlation analy¬ 
sis, discussing the possibilities that its form is attributable 
to selection ehects or dust. In Section 5 we interpret our 
results as indicative of statistical gravitational lensing and 
discuss the implications of such an interpretation. Section 6 
presents our conclusions. 


2 THE 2DF QSO REDSHIFT SURVEY 

The 2dF QSO Redshift Survey (henceforth 2QZ) is named 
for the 2-degree Field multi-object spectrograph being 
utilised at the AAT to survey 740 deg^ of sky. The 2QZ 
patch is contained within the 2dF Galaxy Redshift Survey 
(Colless 1998), and comprises 75° in ascension and 5° in dec¬ 
lination in the region of both the North and South Galactic 
Gaps. In the North, the strip is bounded by —2?5 and 2?5 
Declination, and 9^ 50™ and 14^ 50™ Right Ascension. The 
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southern strip extends from —32?5 to —27?5 Declination and 
from 21^ 40™ to 3^ 15™ Right Ascension. Missing area can 
be attributed to holes containing cuts for bright stars. 

The QSO target sample, the input catalogue (Smith et 
al. 1997), was selected from colour cuts in the (u—hj) against 
(bj — r) plane of APM scanned UKST data to limiting mag¬ 
nitude 18.25 < bj < 20.85 (18.4 ^ b ^ 21). At the bright 
end, the survey has been extended to bj > 16 using the 
6-degree Field spectrograph. Spectra of each object in the 
input catalogue were taken, to determine what percentage 
are genuinely QSOs. 2QZ spectroscopic observations have 
been carried out since 1997, with a projected identification 
of 25,000 QSOs. Each 2dF QSO survey held was observed by 
125 optical hbres for around an hour. Many of the remain¬ 
ing 275 hbres concurrently observed galaxies for the Galaxy 
Redshift Survey. The minimum hbre separation was of the 
order of 20 arcsec on the sky. The spectra cover the wave¬ 
length range 3700-8000A. The process yielded an average 
signal-to-noise ratio of around 5 in the central wavelength 
range of the faintest sources, allowing categorical spectro¬ 
scopic identihcation of objects in all but the poorest seeing. 
For further general information on the technicalities of the 
survey, consider Boyle et al. (2000) or Groom et al. (2001). 

The 2QZ is now essentially complete. Spectra have been 
measured for nearly 45,000 objects, of which about 12,000 
are categorically identihed as stars (or White Dwarfs) and 
about 23,000 as QSOs. The catalogue also contains a number 
of Narrow Emission Line Galaxies (NELGs), BL Lac objects 
and unidentihed, or low signal-to-noise, spectra. The equa¬ 
torial strip near the NGC (North Galactic Gap) contains 
nearly 10,500 of the QSOs. The strip in the SGC (South 
Galactic Gap) contains over 12,800 QSOs. 


2.1 Number counts by magnitude 

The expected strength of lensing-induced correlations be¬ 
tween galaxies and a magnitude-limited sample of QSOs 
depends on the slope of the integrated number-magnitude 
counts, P, fainter than the QSO sample’s limit (Narayan 
1989). We thus need to estimate this slope to interpret our 
results. As we do not have information fainter than the limit 
of the 2QZ, we fit models to the 2QZ number-counts brighter 
than the bj — 20.85 limit and extrapolate the counts to 
fainter magnitudes. Such models are better constrained by 
htting them to the differential counts and integrating them. 

In Figure 1 we present the differential number counts 
by magnitude of the completed 2QZ, which have been ex¬ 
tensively corrected for incompleteness and dust (see Boyle 
et al., in preparation) and averaged over both hemispheres. 
Also plotted are points from the 6dF QSO Survey (6QZ). 
When determining the number-count relation, the sample is 
restricted to the redshift range 0.3 < « < 2.2, the range for 
which the 2QZ is designed to be photometrically complete 
(Groom et al. 2001). Error bars are Poisson. 

The dashed line is a smoothed power law model, where 
the differential counts are expressed in the form 


dN _ No 

dm 10 ““d('"-'"o) -(- 10“/3d(m-mo) 


( 1 ) 


The best-fit model has a bright-end slope of aa = 0.98, 
a knee at mo = 19.1 and a faint-end slope of Pd = 0.15. 
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Figure 1. The differential number counts of the 2dF QSO Red¬ 
shift Survey. The points are QSO number counts in 0.2 mag bins, 
with Poisson errors. The dashed line ia a smoothed power law fit 
to the data. Brighter data points are from the 6dF QSO Redshift 
Survey. Also displayed are the faint data from Boyle, Jones Sz 
Shanks (1991) and Koo & Kron (1988). 
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Figure 2. The N(< m) relationship, or integrated number 
counts, for the 2dF QSO Redshift Survey. The points are QSO 
number counts in 0.2 mag bins, with Poisson errors. The line is a 
smoothed power law fit to the differential number counts. Brighter 
data points are from the 6dF QSO Redshift Survey. Also displayed 
are the faint data from Boyle, Jones Sz Shanks (1991) and Koo Sz 
Kron (1988), which have been offset slightly to prevent the points 
from merging. 
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Figure 3. The N(z) relationship for the 2dF QSO Redshift Sur¬ 
vey. The wider histogram represents QSO number counts in 0.05 
redshift bins. The narrow histogram is an analytic model of the 
galaxy distribution at the B = 20.5 limit of the galaxy samples 
used in this paper. The distributions are normalised to reach a 
maximum at unity. 

This model is consistent with faint data from Boyle, Jones 
& Shanks (1991) and Koo & Kron (1988), which are also 
marked in Figure 1. Data from the literature have undergone 
a zeropoint correction of bj — b — 0.15 and a dust correction 
of about 0.1 mag. 

In Figure 2, the models are integrated and displayed 
against the integrated number-magnitude counts. The best- 
fit model to the differential counts has an average integrated 
faint-end slope of /3 = 0.29. Although the model is well con¬ 
strained by the data (cr = ±0.015) the many incompleteness 
corrections to the faint-end data (again, see Boyle et ah, in 
preparation) mean the la error may be as high as ±0.05. 
When BFS88 modelled the magnitude distribution of QSOs 
as a broken power law, they determined a B-band faint end 
slope of 0.32 — 0.33. In an extensive review, Hartwick & 
Schade (1990) subsequently determined a faint-end slope of 
0.31. Our average slope is thus consistent with these earlier 
authors. 


2.2 Number counts by redshift 

When making model predictions to interpret our results, we 
shall need to assign both QSOs and galaxies a redshift, in 
order to estimate the average QSO-galaxy separation (for a 
given cosmology). We randomly select QSO redshifts from 
the number-redshift distribution of the 2QZ and galaxy red- 
shifts from the analytic distribution of Baugh & Efstathiou 
(1993). The average QSO-galaxy separation is then esti¬ 
mated as the mean separation from thousands of random 
pairs of redshifts. 

In Figure 3 we display the number counts by redshift in 
the 2QZ as a histogram. Also marked is an analytic model 
for a galaxy redshift distribution (Baugh & Efstathiou 1993) 


integrated to S = 20.5, the galaxy sample limit we typically 
consider in this paper. Both distributions have been nor¬ 
malised to peak at 1. Note that for redshifts greater than 
0.4, less than 0.4 per cent of the projected galaxy distribu¬ 
tion overlaps the QSO distribution. 

3 APM AND SDSS CLUSTERS 

3.1 The catalogues 

Our correlation analysis in subsequent sections relies on 
galaxies from two generations of surveys, the APM Galaxy 
Survey and the Early Data Release (EDR) of the Sloan Dig¬ 
ital Sky Survey (SDSS). 

The initial APM survey (Maddox et al. 1990a) was de¬ 
rived from Automated Plate Measuring scans of 185 pho¬ 
tographic UKST plates and covered about 10 per cent of 
the entire sky around the SGP. The original APM region 
was bounded between a Right Ascension of roughly 21*' to 
5^, with Declination from —72° to —18°. Images were de¬ 
tected to a B-band magnitude of 21.5, allowing galaxies to 
be definitively identified down to B < 20.5. The photometry 
was extensively aligned using overlapping plates (Maddox et 
al. 1990b). The APM study was later extended and provided 
the input catalogues for the 2dF surveys. The original APM 
completely covers the southern 2QZ strip. 

The SDSS (http://www.sdss.org) is imaging the north¬ 
ern sky in five bands designed for GCD photometry 
{u',g',r',i' and z'). The survey will trace an ellipse cen¬ 
tred on 12** 20™ Right Ascension and 32.8° Declination, 
roughly extending from 7^ 6™ to 17^ 34™ in Right Ascension 
and ±55° of declination. The SDSS should be complete to 
g' ~ 23.3 and r' ~ 23.1 (York et al. 2000), about equivalent 
to 6 ~ 23.5 using a typical colour transformation (Yasuda et 
al. 2001). The Early Data Release (Stoughton et al. 2002) 
contains about the first 460 deg^ of the SDSS. Conveniently, 
one strip, extending from 9*“ 40™ to 15*“ 44™ in Right As¬ 
cension and —1° to 1.5° Declination, is contained within the 
northern 2QZ strip. 

3.2 Objective group catalogues 

We wish to correlate the 2QZ sample with groups of galaxies 
rather than the field as any correlation attributed to leasing 
will be stronger for more massive foreground structures. A 
catalogue of objectively identified groups in the Southern 
APM is available in the literature (Dalton et al. 1997) but 
comparatively little has been available in the region of the 
northern 2QZ since the publication of the AGO catalogue 
(Abell, Corwin & Olowin 1989). Additionally, it is useful 
to have listings of the position of each galaxy within each 
cluster, not merely the positions of cluster centres. Eollowing 
BFS88 we turn to the clustering routine of Turner & Gott 
(1976) to objectively identify clusters of galaxies in the APM 
survey and SDSS EDR. 

To determine the boundaries for cluster membership 
in a galaxy sample, we assign an overdensity, S, a factor 
by which we wish our group density to exceed the mean 
surface density across the entire region we consider (d). We 
then calculate the largest possible circle of radius 9c such 
that 
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Figure 4. Objectively defined galaxy groups in the Sloan Digital Sky Survey Early Data Release. Each point is a galaxy in a region 
sufficiently dense that it meets clustering criteria outlined in the text. The axes correspond to the limits of the 2dF QSO Redshift Survey. 
The coordinate system has been transformed to B1950. 


a{6 < 9c) > 5(7 (2) 

where cr is the surface density of galaxies within the (cir¬ 
cular) region centred on any particnlar galaxy in our sample 
and enclosed by an angular radius 9. Over the miniscule 
angles typically considered, a can be expressed as 

^ ^ N{< 9) ^ N{< 9) 

27r(l —COS0) n9^ 

where N{< 9) is the total number of galaxies within an 
angle 9 of the particular galaxy we are considering, including 
the particular galaxy itself. The critical angular radius 9c 
is determined for each individual galaxy in the sample. A 
group is defined as all galaxies that have overlapping critical 
radii. 

Two classification questions remain; how overdense are 
groups (what value should 5 take)? and what size must 
groups attain before we call them a group (what is the min¬ 
imum number of galaxies, Nmin, in a group)? We take the 
values of (5 = 8 and Nmin = 7 chosen by Stevenson, Fong 
& Shanks (1988) in a similar analysis, and used by BFS88. 


The choice of the overdensity parameter was originally sug¬ 
gested by Turner & Gott (1976) and weighs the possibility 
of losing poor clusters at high 5 values against false grouping 
of galaxies in the field at low 5. The choice of A = 7 as a 
minimum group size reduces the likelihood of chance align¬ 
ments of galaxies at different redshifts being grouped. The 
parameters reflect a distribution of groups that may have 
been selected by eye (Stevenson, Fong & Shanks 1988). 

Figures 4 and 5 show the groups determined from the 
SDSS EDR data in the northern strip and from APM data 
in the southern strip. In both cases, the galaxy sample was 
limited to B < 20.5, the limit of the APM. The SDSS EDR 
data was transformed from the SDSS photometric system 
using the colour transformation of Yasuda et al. (2001). The 
SDSS EDR coordinate system has been transformed from 
J2000 to B1950, the system we shall use throughout this 
paper. The axes in both figures correspond to the limits of 
the 2QZ. A great deal of the structure produced by Dal¬ 
ton et al. (1997) is reproduced well in our southern group 
catalogue, even though Dalton used different clustering cri¬ 
teria. 30 per cent of Dalton’s cluster centres within the 2QZ 
boundary lie within 2 arcmin of our group centres, improv¬ 
ing to 75 per cent of Dalton’s cluster centres lying within 
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Figure 5. Objectively defined galaxy groups in the Southern APM Survey. Each point is a galaxy in a region sufficiently dense that it 
meets clustering criteria outlined in the text. The axes correspond to the limits of the 2dF QSO Redshift Survey. The coordinate system 
is B1950. 


5 arcmin of our own. 90 per cent of group centres match up 
within a separation of 7.7 arcmin, the mean cluster radius 
of Dalton’s groups within the 2QZ SGC region. Many of 
the larger groups automatically picked out by our technique 
also correspond to previously identified rich clusters (Abell 
1958), for instance, the large cluster around 14^ 12™ Right 
Ascension in Figure 4 corresponds to Abell 1882 (Richness 
Class III, 2 ~ 0.137). Abell assigned 166 galaxies to Abell 
1882, we determine 153 members. Our voids are also in good 
agreement with the literature. The sparse region we pick out 
in the Southern APM from Right Ascension 2*“ 30™ to 3^ 
appears to be real. In this region, Dalton (1997) finds no 
clusters over the Declination range —32?5 to —28° and Abell 
finds a single (Richness Class 0) cluster in the range —32?5 
to —29? 5. The void is not an effect of intervening dust ob¬ 
scuring galaxies out of the APM - the region of the SCC 
strip from 23*“ to 23^ 30™ has more dust across it (accord¬ 
ing to the maps of Schlegel, Finkbeiner & Davis 1998) and 
a denser population of groups. In any case, we note that 
similar voids are picked out at, say, 22^ 30™ in the Southern 
APM, or 10'’ 30™ in the SDSS EDR. 


4 CROSS-CORRELATION ANALYSIS 

We now turn our attention to cross-correlations of galaxy 
groups with 2QZ QSOs. 


4.1 Method 

The two-point angular correlation function uj{0) measures 
the probability (dP) of finding pairs of sources of mean num¬ 
ber density n within solid angle dQ separated by angle 9 

dP = hih2(l-I-cu(6'))df2idl22 (4) 

The correlation function then, measures the relative 
amount of associated structure in pairs of sources for var¬ 
ious angular separations. If there is some association the 
correlation function will be positive. If there is some avoid¬ 
ance, the correlation function will be negative. If there is no 
distributed pattern whatsoever, as we might expect from, 
say, stars versus extragalactic objects, the correlation func¬ 
tion will be zero. Throughout this analysis, we estimate the 
two-point correlation function using an equation of the form 
proposed by Peebles (1980) 
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uj{e) 


DD{6)nB. 

DR(6)nr> 


(5) 


where D refers to a data point (either galaxy, QSO 
or star) and R refers to a random point (mock galaxy or 
QSO). DD, then, could be the total number of galaxy- 
QSO, QSO-galaxy, galaxy-star or star-galaxy pairs, to name 
a few constructions. If DD were a reference to the number of 
QSO-galaxy pairs, then DR would be the number of QSO- 
random pairs, where the random catalogue consisted of a 
mock galaxy distribution. The parameters fin and fin re¬ 
fer to the total number density of data points and the total 
number density of points derived from a random catalogue, 
respectively. The term 6 refers to the angular separation be¬ 
tween the pairs. 

Random catalogues are constructed by randomly creat¬ 
ing points with the same overall distribution as the sources. 
For example, a random galaxy catalogue contains points 
randomly distributed across the sky except where there are 
holes in the ARM or SDSS and a random QSO catalogue 
contains points distributed according to the coverage map 
of the 2QZ. The 2QZ coverage map is produced by calcu¬ 
lating the ratio of spectroscopically observed objects to in¬ 
put catalogue objects in each region defined by the inter¬ 
section of 2dF pointings. This distribution is then pixelized 
using arcmin^ pixels (Groom et al. 2001). When estimating 
the correlation function, the random catalogue is always at 
least 50 times bigger than the data catalogue. Note that the 
cross-correlation between, say, galaxies and QSOs, could be 
performed in one of two directions, by centring on QSOs 
and counting galaxies or by centring on galaxies and count¬ 
ing QSOs. With ideal samples, these procedures should be 
equivalent. In many cases, if gradients or biases that are not 
accounted for in the random catalogue exist in one or both 
samples, the two directions may not be exactly equivalent. 

Since we are mainly interested in small scale (~ 10 ar- 
cmin) cross-correlations, we measure them locally by split¬ 
ting both the NGC and SGG into 15 individual fields. In 
the north, this is purely arbitrary but in the south, the field 
boundaries correspond to the edges of the APM plates. Gor- 
relations are then counted within each individual held, with 
the resulting numbers of pairs being totalled for all helds to 
yield the total number of pairs across the entire survey area. 
The total number of data-data (DD) pairs and data-random 
(DR) pairs are then taken in ratio, as per Equation (5), to 
estimate the global correlation function. This held-to-held 
analysis should nullify the effects of different photometric 
zero-point calibrations between plates, or gross variations 
across strips. Errors are determined via the spread in the lo¬ 
cal correlation function between each held. Essentially, these 
held-to-held errors are Icr standard deviations in the value 
of the correlation function between helds, inverse variance- 
weighted to account for the different numbers of sources on 
each held. Errors are thus estimated via 
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where is the held-to-held error on the correlation 
function, and the subscript L stands for local, referring to 
an individual plate. The DRl/DR factor, then, weights each 
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Figure 6. The cross-correlation between 2QZ QSOs, and galax¬ 
ies of limiting magnitude B = 20.5 found in objectively derived 
groups of at least seven members, for both 2QZ strips. Crosses 
correspond to the NGC strip, triangles to the SGC strip. The 
numbers of each sample within the boundaries of the 2QZ are 
displayed. Errors are field-to-field. 




Figure 7. The cross-correlation between galaxies of limiting mag¬ 
nitude B = 20.5 found in objectively derived groups of at least 
seven members and 2QZ QSOs, combined for both the south¬ 
ern and northern 2QZ strips. Both results centring on QSOs and 
counting galaxies (triangles) and centring on galaxies and count¬ 
ing QSOs (crosses) are presented. The numbers of each sample 
within the boundaries of the 2QZ are displayed. Errors are field- 
to-field. 
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Figure 8. The comparison between Figure 7 and the result of 
BFS88. The BFS88 data (open circles) have been scaled to ac¬ 
count for their stellar contamination. 


field so that fields with more objects are lent more signifi¬ 
cance in the error calcnlation. 

We have compared this field-to-field error estimate to 
Poisson error, as estimated by (1 -|- uj{9))/-\/DD. Although 
the two types of error estimate show a strong relation, the 
Poisson estimate overestimates the error, as compared to 
the variance in random simulations of the data. This might 
be expected as an individual source may be associated with 
structure at several separations, so errors over a given range 
can be strongly correlated. On the smallest scales, especially 
on scales where there are no data-data pairs on some of 
the plates, even the field-to-field estimate of error breaks 
down. When fitting models to the data, we use an estimate 
of the error based on the variance between many random 
simulations of the data weighted to reproduce the model of 
interest. 

In the following subsections, we calculate a variety of 
field-to-field cross-correlation functions using Equation (5) 
for 2QZ QSOs against groups of APM or SDSS galaxies. 

4.2 Cross-correlation of QSOs and group galaxies 

In Figures 6 (NGC and SGC) and 7 (combined) we dis¬ 
play the cross-correlations between spectroscopically iden- 
tifed 2QZ QSOs and galaxies in groups of at least 7 mem¬ 
bers objectively derived from the SDSS EDR in the NGG 
strip and the APM catalogue, in the SGG strip. In Figure 7, 
we show both directions of correlation, we have centred on 
galaxies in groups and counted QSOs, and have centred on 
QSOs and counted galaxies. Figure 8 shows the comparison 
between our result and BFS88, where the BFS88 data have 
been scaled to allow for their projected 25 per cent contam¬ 
ination by stars. The results are displayed at the smallest 
scales for which we still believe the field-to-field-errors. The 
largest scale displayed is a few bins before edge effects caused 
by the 2QZ strips only being 300 arcmin in declination begin 



Figure 9. The cross-correlation between 2QZ stars, and galax¬ 
ies of limiting magnitude B = 20.5 found in objectively derived 
groups of at least seven members, for both 2QZ strips. Crosses 
correspond to the NGC strip, triangles to the SGC strip. The 
numbers of each sample within the boundaries of the 2QZ are 
displayed. Errors are field-to-field. 


to have any effect. The numbers displayed on each plot are 
the numbers of objects from each sample present within the 
confines of the 2QZ boundary. Note that the SDSS EDR 
sample is half the size of the northern 2QZ strip, so does 
not contribute as significantly to the combined signal as the 
APM. A redshift cut has been made at z = 0.4 in the QSO 
sample; if the signal is caused by lensing, this will theo¬ 
retically reduce the overlap between QSOs and foreground 
matter to at most 0.4 per cent (see Section 2, above). 

Figure 6 shows that there is good consistency between 
the northern and southern correlation functions for QSOs 
versus galaxies in groups, justifying combining the signals. 
Figure 7 demonstrates that there is excellent consistency in 
the cross-correlation signal between QSOs and group galax¬ 
ies irrespective of the direction in which the function is cal¬ 
culated, suggesting the signal is robust, free from the in¬ 
fluence of any gradient or incompleteness in the samples 
used. There is a significant 3 .Oct anti-correlation between 
galaxy groups and spectroscopically identified 2QZ QSOs 
on scales out to 10 arcmin, based on collecting the data in 
a single 10 arcmin bin and calculating the rms field-to-field 
variation. For an Einstein-de-Sitter cosmology, 10 arcmin 
would translate to 1 h~^ Mpc at an average group redshift 
of « ~ 0.15 (1.1 h~^ Mpc for ACDM. The average anti¬ 
correlation in such a 10 arcmin bin is —0.049. Our data 
compare well in Figure 8 with the anti-correlation discov¬ 
ered by BFS88 when correlating a UVX object sample with 
cluster galaxies. BFS88 declared a more significant 4 ct sig¬ 
nal on <4 arcmin scales. We have also measured the same 
effect after correcting the random QSO catalogue for Galac¬ 
tic dust, using Schlegel maps (Schlegel, Finkbeiner & Davis 
1998). Although the strength of the effect is not altered, the 
variation between fields is reduced, meaning the significanee 
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Figure 10. The cross-correlation between galaxies of limiting 
magnitude B = 20.5 found in objectively derived groups of at 
least seven members and 2QZ stars, combined for both the south¬ 
ern and northern 2QZ strips. Both results centring on stars and 
counting galaxies (triangles) and centring on galaxies and count¬ 
ing stars (crosses) are presented. The numbers of each sample 
within the boundaries of the 2QZ are displayed. Errors are field- 
to-field. 

of the result increases slightly to 3.3a. Note that, as neither 
of the galaxy populations we use are corrected for Galac¬ 
tic dust, no correction for Galactic dust in the QSO sample 
should be necessary when cross-correlating the 2QZ with the 
galaxy populations. In any case, correcting for Galactic dust 
has a very small effect. 

To test the anti-correlation between group galaxies and 
QSOs, we need a control sample that has been through 
the same measurement and reduction as the QSO sample. 
BFS88 found no correlation between a large (~ 27,000) con¬ 
trol sample of non-UVX stars and group galaxies, however, 
our situation is notably more complicated. As we require a 
control sample that has been through the same processes 
as the QSOs, including spectroscopic identification, we are 
restricted to a smaller sample of stars (~ 12,000) compared 
to our QSO sample (~ 22,500). Additionally, the selection 
criterion of the 2QZ tends to pick out specific populations 
of stars, such as White Dwarfs. 

In Figures 9 (NGC and SGC) and 10 (combined) we 
show the cross-correlation between spectroscopically identi¬ 
fied 2QZ stars and galaxies in groups. Unlike the QSO result, 
the NGC and SGC cross-correlations are not entirely con¬ 
sistent, suggesting that the different physical distribution of 
stars in the two strips effects the correlation function. The 
combined results are also inconsistent, although the errors 
are sizable. The combined results show evidence of gradi¬ 
ents on large scales, as, unlike in the case of the QSOs, the 
cross-correlation result is dependent on whether we centre 
on stars and count galaxies or centre on galaxies and count 
stars. Indeed, the distribution of stars within our galaxy does 
display a gradient with galactic latitude in both the NGC 
and SGC, unlike the distribution of QSOs, which is flat. 


The cross-correlation between stars and galaxies has some 
negative points on 3-7 arcmin scales but does not have the 
same form as the QSO anti-correlation. In particular, the 
star-signal does not continue to decrease on scales less than 
3 arcmin. If we calculate the significance of the signal out 
to 10 arcmin, we find a 1.6a (1.4a with correction for dust) 
anti-correlation for the result centring on stars and count¬ 
ing galaxies and a 3.0a (2.9a with dust-correction) anti¬ 
correlation for the result centring on galaxies and counting 
stars. Much of this discrepancy is caused by the large-scale 
gradients; if the galaxy-centred result is moved up so there 
is no anti-correlation on large scales, it comes into line with 
the star-centred result. The anti-correlation in both of the 
combined star results is caused entirely by the few points on 
3-7 arcmin scales. Surprisingly, most of the anti-correlation 
between stars and group galaxies is actually caused by White 
Dwarfs. If we discard the White Dwarfs from the stellar sam¬ 
ple, there is no significant anti-correlation between stars and 
group galaxies (see Figure 11). Whether there is a physical 
reason for this (circumstellar dust?) remains to be seen, oth¬ 
erwise the most likely explanation for the anti-correlation 
seen at 3-7 arcmin in both the stars and White Dwarfs is a 
statistical fluctuation due to the low numbers of these ob¬ 
jects. 

We regard the results using the above control samples 
of stars as encouraging in terms of ruling out a systematic 
source for the QSO-group anti-correlation. However, because 
of the low numbers of stars, their gradients in Right Ascen¬ 
sion and the anti-correlation detected at 3-7 arcmin, stars 
may not form the ideal control sample and there may be a 
residual doubt as to whether there is a systematic contribu¬ 
tion to the QSO anti-correlation on 3-7 arcmin scales caused 
by the fibre positioning constraint of the 2dF instrument. 
After all, 2dF observed the b < 19.5 galaxies and QSOs si¬ 
multaneously and in dense fields close pairs of objects may 
have been missed due to the minimum fibre spacing, even 
though 2dF candidates were given a higher priority in the 
fibre allocation to prevent imprinting the galaxy structure 
on the QSO distribution. It is easy to show there is no fibre 
positioning effect by comparing the cross-correlation of the 
46,000 objects in the 2dF input catalogue with the 43,000 
objects that were observed spectroscopically. These results 
are shown respectively in the upper left-hand and lower left- 
hand panels of Figure 11. These two results are in all respects 
identical with no systematic difference between them, lead¬ 
ing us to conclude that there is no anti-correlation induced 
on these scales by the hbre positioning constraint. Therefore 
we conclude that the QSO-group galaxy anti-correlation is 
probably real. 

Note that the reason the correlation between group 
galaxies and the input catalogue is flat, even though it con¬ 
tains the significantly anti-correlated QSO signal, becomes 
apparent when we split the input catalogue up into its con¬ 
stituent parts. The right-hand four panels in Figure 11 dis¬ 
play the main subsamples of the input catalogue. We can see 
that the Narrow Emission Line Galaxies included in the in¬ 
put catalogue exhibit a positive correlation that cancels out 
the anti-correlation exhibited by QSOs, leaving the input 
catalogue result flat. 
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Figure 11. The cross-correlation between 20204 galaxies of limiting magnitude B = 20.5 found in objectively derived groups of at least 
seven members and samples of objects taken from the 2QZ, combined for both the northern and southern 2QZ strips. The upper-left 
panel is the input catalogue. The lower-left panel is the fibre catalogue. The other panels are specific objects that comprise the majority 
of the fibre catalogue; stars, QSOs, White Dwarfs and NELGs. The numbers of each subsample are displayed. 


4.3 The dust hypothesis 

BFS88 originally attributed the anti-correlation between 
QSOs and galaxies in groups to interim dust in clusters, 
finding that an absorption in the 6—band of Ah ~ 0.2 mag 
was sufficient to cause their observed anti-correlation. Such 
absorption is at the upper limit allowed by Ferguson (1993), 
who quoted a maximum reddening of E{B — V)< 0.06 from 
a composite study of the Mg 2 index of 19 nearby clusters 
and rich groups. Using colour information provided by the 
2dF survey, we can limit the culpability of dust in caus¬ 
ing the anti-correlation signal. The 2QZ measures u — bj 
and bj — r colours. If the anti-correlation between QSOs and 
group galaxies were due to dust, we would expect to observe 
a complementary reddening of QSOs. 

In Figures 12 and 13 we show the distributions of 
colours of 2QZ QSOs that lie within 10 arcmin of any group 
centre (‘obscured’ QSOs) and QSOs that do not lie within 
10 arcmin of any group centre (‘unobscured’ QSOs). Fig¬ 
ure 12 is for u — bj colours. Figure 13 for bj — r colours. 
We have selected 10 arcmin as the radius of interest as it 


corresponds to the extent of the anti-correlation signal in 
Figure 7. Were the anti-correlation due to intervening dust 
in galaxy groups, as proposed by BFS88, we would expect 
to see a complementary reddening of QSOs on the scale of 
the anti-correlation, and the distributions of ‘obscured’ and 
‘unobscured’ QSOs would differ. 

A two-sample Mann-Whitney U-test fails to reject the 
null hypothesis that the ‘obscured’ and ‘unobscured’ dis¬ 
tributions have the same median, for both the {bj — r) and 
(u—bj) colours. A two-sample Kolmogorov-Smirnov test fails 
to reject the null hypothesis that the ‘obscured’ and ‘unob¬ 
scured’ distributions are drawn from the same parent popu¬ 
lation, again for both colours. Student’s t-test demonstrates 
that the means of the ‘obscured’ and ‘unobscured’ samples 
are in agreement. For the (u — bj) colour distributions, the 
‘unobscured’ mean and standard error are —0.6789 ± 0.0030 
and for the ‘obscured’ are —0.6687± 0.0064. For the (bj — r) 
colour distributions, the ‘unobscured ’ mean and standard 
error are 0.3644 ±0.0030 and the ‘obscured’ mean and stan¬ 
dard error are 0.3626 ± 0.0065. There are 4025 ‘obscured’ 
QSOs and 17752 ‘unobscured QSOs. Student’s t-distribution 
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Figure 12. The relative distribution of u — bj colours of 2QZ 
QSOs. The solid histogram shows the colours of QSOs that do 
not lie within 10 arcmin of a group centre. The dashed histogram 
shows the colours of QSOs that lie within 10 arcmin of a group 
centre. 
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Figure 13. The relative distribution of — r colours of 2QZ 
QSOs. The solid histogram shows the colours of QSOs that do 
not lie within 10 arcmin of a group galaxy. The dashed histogram 
shows the colours of QSOs that lie within 10 arcmin of a group 
galaxy. 



sets the following 95 (99) per cent upper limits on red¬ 
dening between the two distributions: E{u — bj) = 0.012 
(E(u - bj) = 0.016), E(bj - r) = 0.012 {E{hj - r) = 0.016). 
Now, our average group size is around 2.5 arcmin and larger 
groups in our sample have an angular size of 5 arcmin, so we 
might also be interested in any reddening on these scales. 
Repeating the above analysis on these scales. Student’s t- 


test suggests the 95 (99) per cent reddening limits between 
the QSO population within 5 arcmin of any group centre 
and the QSO population lying greater than 10 arcmin are 
E{bj — r) = 0.020 (0.028) and the 95 (99) per cent redden¬ 
ing limits between the QSO population within 2.5 arcmin of 
any group centre and the remainder are E{bj — r) = 0.039 
(0.056). These reddening limits apply for both the u — bj and 
bj—r colours. The limits inevitably increase as the ‘obscured’ 
population dwindles in size, although there are still 286 
QSOs within 2.5 arcmin of any group centre. Again, for both 
the 2.5 arcmin and 5 arcmin scales, the Kolmogorov-Smirnov 
test failed to reject the hypothesis that the ‘obscured’ and 
‘unobscured’ distributions are drawn from the same parent 
distribution of colours. Assuming the usual Galactic dust 
law, the 95 per cent limit from the bj—r colours on B-band 
absorption within 5 arcmin (~ 0.5 h~^ Mpc at 2 ~ 0.15) of 
any group centre is Ab < 0.06 mag. The similar limit within 
2.5 arcmin of any group centre is Ab < 0.13 mag. A simple 
model (see BFS88 for details), taking the slope in the QSO 
number magnitude counts to be 0.29, suggests that these 
levels of absorption correspond to a correlation function of 
oJcq < —0.039 for As < 0.06 mag and Ucq < —0.077 for 
As < 0.13 mag. For the entire sample out to 10 arcmin, the 
95 per cent reddening limits suggest an anti-correlation of 
only —0.029. 

It seems that 2QZ QSOs close to galaxy groups are in¬ 
sufficiently reddened to explain the anti-correlation signal as 
an effect of dust in galaxy groups. It might be argued that as 
the 2QZ is colour-selected, our reddening values are biased 
by the limits set on the colours of 2QZ QSOs. We do not 
mean to argue that our values are objective determinations 
of the extent of dust in galaxy groups, only that there is in¬ 
sufficient reddening of QSOs close to galaxy groups within 
the 2QZ to explain the anti-correlation signal. Undoubtedly, 
there will be some heavily reddened QSOs close to group 
centres that the 2QZ fails to observe, however, the QSOs 
the survey does observe have no tendency to redness close 
to group galaxies. The 2QZ sample colour-colour distribu¬ 
tion peaks significantly bluewards of any colour limit, so we 
believe that the low reddening measure is not forced by the 
colour selection procedure. 

Previous evidence for dust in galaxy groups has been 
controversial, with some authors claiming detections and 
others claiming upper limits. Girardi et al. (1992) confirm 
galaxies in local groups that are blueshifted relative to the 
group average tend to have a larger colour excess than that 
group average, and suggest background galaxies falling to¬ 
wards the group centre suffer reddening by dust in the group. 
The amount of dust Giardi et al. (1992) infer is E{B — V) 
~ 0.1 — 0.2 on 0.75 Mpc scales, equivalent to an absorp¬ 
tion of Ab ~ 0.4 — 0.8 mag (for the usual galactic dust 
law). Such an amount of absorption is at odds with the 
upper limit quoted by Ferguson (1993). Ferguson studied 
the correlation of Mg2 index with B — V colour for local 
groups and clusters, and compared these results to a sample 
of field galaxies. Ferguson (1993) quoted E{B — U) < 0.06 
as a conservative (90 per cent) upper limit on reddening 
in a sample of clusters and groups and, similarly, an upper 
limit of E{B — U) < 0.05 (Ab < 0.2 mag) for a sample of 
poorer groups. Ferguson (1993) considers a scale of 0.5 Mpc 
to be the central contribution to his reddening limits. The 
95 per cent upper limits on our reddening results out to 
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0.5 h~^ Mpc (5 arcmin) allow for us much absorption by 
dust as Ab < 0.06 mag, more consistent with the results of 
Ferguson (1993) rather than those of Girardi et al. (1992). 

The possibility remains that some specifically tailored 
dust model could still explain the anti-correlation. One pos¬ 
sibility is a smooth distribution of dust around galaxy groups 
that does not obey the reddening laws observed in our 
galaxy (so-called ‘grey dust’). A less ad hoc explanation is 
a patchy distribution of dust around galaxy groups that is 
heavily concentrated in the line of sight direction, such that 
QSOs are completely obscured without being reddened. Av¬ 
eraged over many groups, this would have the effect of re¬ 
moving QSOs from the sample near galaxy groups with¬ 
out overtly reddening the QSO sample on similar scales. Fi¬ 
nally, we must consider the possibility that multiple effects 
of dust and lensing combine to produce the observed anti¬ 
correlation. So, although our current analysis can rule out 
a smooth distribution of typical dust around galaxy groups, 
we cannot deny more tailored dust models, such as grey 
dust, patchy dust or a combination of gravitational lensing 
and dust. 


5 WEAK LENSING 

We have seen that the 3 .Oct anti-correlation between 2QZ 
QSOs and galaxies in groups is unlikely to be either a se¬ 
lection effect or a product of normal, smoothly distributed 
dust in galaxy groups obscuring QSOs. In this section we 
consider the possibility that the anti-correlation results from 
the statistical lensing of QSOs by foreground galaxy groups 
(Groom & Shanks 1999). 

5.1 Modelling approach 

In the Appendix, we outline models we use to describe how 
the statistical lensing of background QSOs may trigger an 
anti-correlation signal of the type discovered above. The 
models describe a slightly different situation than hitherto 
discussed, the lensing of background flux around the centres 
of dark matter profiles, and so describe the correlation func¬ 
tion of QSOs against the centres of galaxy groups. Previ¬ 
ously, to compare our results with BFS88, we have measured 
the correlation between QSOs and individual galaxies in 
groups. To facilitate comparison with models, we now cross¬ 
correlate QSO positions with group centres, rather than 
group galaxies. Additionally, in prior Sections larger groups 
of galaxies have been weighted more highly as each indi¬ 
vidual galaxy position within the group would be counted. 
We fit models to data where the correlation function is un¬ 
weighted - each cluster is considered equally. 

The cross-correlation between 2QZ QSOs and the cen¬ 
tres of galaxy groups objectively determined from the 
APM South and the SDSS EDR is displayed in Figure 14. 
For comparison, both the weighted and unweighted cross¬ 
correlations between QSOs and the centres of galaxy groups 
are shown. The weighted result (crosses) is a reflection of 
the analysis in Section 4.2, the QSOs are counted against 
each galaxy in a group, so larger groups are lent more sig¬ 
nificance. The unweighted result is not biased by the size of 
the groups. Binning the data in a sole bin out to 10 arcmin 
and measuring the rms field-to-field variation in this bin, the 



9 (arcmin) 

Figure 14. The cross-correlation between 2QZ QSOs and ob¬ 
jectively derived galaxy groups of limiting magnitude B = 20.5 
combined for both the southern and northern 2QZ strips. The 
numbers of each sample within the boundaries of the 2QZ are 
displayed. Crosses are derived by weighting each galaxy group by 
the number of galaxies in the group and have been offset slightly 
for display purposes. Triangles are derived counting each group 
equally. Errors are field-to-field. Also marked as a dashed line 
is the unweighted result for only those groups with 14 or fewer 
members. Errors on this line are as plotted for the triangles. 

average anti-correlation of the weighted result is —0.049 with 
a significance of 2.9 ct. Comparing with the results in Section 
4.2 , the strength and significance of the result weighted by 
galaxy number proves essentially identical whether we cor¬ 
relate QSOs with galaxies in groups or the group centres. 
The unweighted result has an average anti-correlation out to 
10 arcmins of —0.034 with a significance of 2.9 ct. When we 
do not weight the cross correlation by the number of galax¬ 
ies in the group, the strength of the result thus drops by 
30 per cent compared to the anti-correlation between QSOs 
and groups outlined in Section 4.2. We can deduce that the 
anti-correlation signal is stronger for larger groups, as would 
be expected if it is due to lensing. In fact, it is straightfor¬ 
ward to calculate that if we assume a linear relationship 
between the mass of a group and the number of galaxies in 
a group we would expect just such a 30 per cent drop. To 
further demonstrate that the cross-correlation signal is not 
biased up by the handful of large clusters in our sample, 
we have cross-correlated QSOs against a sample of galaxy 
groups that has the largest 10 per cent of objects removed. 
The resulting sample of galaxy groups has between 7 and 
14 galaxies per group. The anti-correlation for this sample 
is plotted as the dashed line in Figure 14 and has essentially 
identical error bars to the full (unweighted) sample (marked 
with triangles). There is practically no difference in the re¬ 
sult if we remove the largest 10 per cent of objects from our 
group sample, the anti-correlation signal thus derived has a 
strength of —0.032 and a significance of 2.5 ct. 

We model the lensing groups as dark matter haloes, 
either Singular Isothermal Spheres or NFW profiles (as de- 
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scribed in the Appendix). In the case of the SIS, the free 
parameter is the velocity dispersion of the sphere. The free 
parameter in the NFW model is the mass within Mpc 

of the centre of the halo centre. For the lensing analysis, we 
nse a faint-end QSO number-count slope of (3 — 0.29 (see 
Section 2). The overwhelming majority of QSOs in the 2QZ 
lie fainter than the knee of a broken power-law model (only 
10 per cent are brighter than bj = 19.1. We have repro¬ 
duced models both approximating the number-magnitude 
counts as a single power-law with slope /? = 0.29 and using 
the full smoothed power law determined in Section 2 (see 
Equation 1) and find no significant difference between the 
two approaches. 

The lens and source distances required by the lens¬ 
ing models (see Equation A 6 , for example) are calculated 
by randomly sampling redshifts from the distributions dis¬ 
played in Figure 3 and deducing the average separation of a 
galaxy-QSO pair, as outlined in Section 2. The redshifts are 
translated to distances using an Einstein-de Sitter cosmol¬ 
ogy- 

When fitting models, we use an estimate of the error 
based on 500 ‘mock’ QSO catalogues of the same size and 
completeness as the 2QZ, calculating the cross-correlation 
between these mock catalogues and our galaxy group cen¬ 
tres in the normal way, and then determining the rms er¬ 
rors between the catalogues. We fit a reasonable model to 
the anti-correlation result, then distribute the random place¬ 
ment of QSOs in the mock catalogues to reflect that model. 
So, if the cross-correlation fit has a value of -0.1 at 2 arcmin, 
a random QSO that is generated 2 arcmin from a group cen¬ 
tre is only 90 per cent as likely to be included in the mock 
QSO catalogue as one a long way from any group. Similarly, 
a QSO that then lies within 2 arcmin of 2 group centres is 
81 per cent as likely to be included. We have tested the in¬ 
dependence of these errors by measuring the covariance of 
points in adjacent bins averaged over the 500 simulations, 
finding the covariance insignificant on all scales. 

5.2 Model fits 

Figure 15 shows the best-fit models for the SIS and NFW 
lensing haloes obtained by minimising the statistic. The 
errors are calculated from the standard deviation in the anti¬ 
correlation signal of 500 mock QSO catalogues as outlined 
above, with ‘mock’ QSOs distributed according to a rea¬ 
sonable model. In the SIS case, mock QSOs are distributed 
according to a a = 600 km s“^ model when calculating er¬ 
rors. In the NFW case, mock QSOs are distributed according 
to a Ml,5 = 3 X 10^^ h~^ Mq model when calculating er¬ 
rors. For brevity, the models are both displayed in Figure 
15 against the error bars calculated for the NFW case. We 
determine the best-fitting models out to 10 arcmin, the ex¬ 
tent of the anti-correlation. Once the best-fitting model is 
determined, the errors are scaled so the reduced is equal 
to 1 and then Icr error-bars on the best-fitting model are 
calculated from this renormalised distribution. The best- 
fit SIS has a velocity dispersion of cr = 1156 ±327 km s“^ 
with a reduced of 0.8. The best-fit NFW has a mass of 
Ml,5 = 1.2 ± 0.9 X 10^® h~^ M0 with a reduced x^ of 1.2. 
The data cannot distinguish variations in the NFW 7 pa¬ 
rameter in the range | < 7 < 1 (see Equation All). As 
7 increases, the predicted anti-correlation decreases below 


1 arcmin and a test between these cases shonld be possible 
in bigger datasets. In the current datasets there is hardly 
enough power to distinguish between the best-fitting SIS 
and NFW models. 

The Einstein Radius of the best-fitting SIS model is 
aronnd 30 arcsec, which results in a radical dip in the solid 
line in Figure 15, corresponding to the terms in the denom¬ 
inator of Equation A 6 being equal. At separations within 
the Einstein Radius, we enter the strong lensing regime and 
the SIS model predicts that each source QSO will produce 
two images. One of these images has already been covered 
in the SIS model prediction for scales larger than the Ein¬ 
stein Radius. The second image appears within the Einstein 
Radius and might be considered a further prediction of the 
SIS model. The divergent nature of the SIS density profile 
on small scales means the dip in Figure 15 may be predicted 
to occur at a scale which is unphysically large. Details of 
the model near the Einstein Radius do not affect the fitted 
SIS masses much, as evidenced by their similarity to those 
masses determined from NFW fitting where the Einstein Ra¬ 
dius is much smaller. The full consequences of strong lensing 
in the 2QZ, including an analysis of the numbers of multiply- 
lensed QSO systems, are discussed elsewhere (by Miller et 
al. 2003). 

As most of our model analysis is made in the weak 
lensing regime, we might be wary of any fit to the small¬ 
est scale data points in the SIS case. As a consistency check, 
we have used mock QSO catalogues to make a direct test 
of the significance of rejection of a = 600 km s”*^ and 
a — 1140 km s“^ SIS models, since the mock catalogues 
were produced for these specific cases. We determine how 
often cross-correlating mock QSO catalogues with galaxy 
groups can return an anti-correlation of significance —2.9 ct, 
as found for the real 2QZ data. We display this result in 
Figure 16. We have created 250 mock QSO catalogues with 
the same size and completeness as the 2QZ. These are then 
cross-correlated against our galaxy groups and the strength 
of the cross-correlation is measured for each mock QSO cat¬ 
alogue. The measure of the significance of a cross-correlation 
is as we have used throughout this paper, based on counting 
the data in one bin out to 10 arcmin. The mock QSO cata¬ 
logues may also be distributed, as outlined above, to reflect 
various models for the lensing halo. We find that an anti¬ 
correlation of —2.9 ct is measured 0/250 times if there are no 
lensing haloes, only 7/250 times (2.2 ct) if the lensing haloes 
are modelled as an SIS with a 600 km s“^ velocity disper¬ 
sion (roughly equivalent to Qm = 0.3) and 137/250 times 
for an SIS with an 1140 km s“^ velocity dispersion (roughly 
equivalent to Qrn = 1). If we had chosen to plot the strength 
of the anti-correlation rather than the significance, we would 
find that a 600 km s“^ model produces an anti-correlation 
of strength -0.034 (as found for the data) only 13/250 times 
( 1 . 9 ( 7 ). Either way, these results reject the 600 km s“^ model 
at about the 5 per cent significance level. Note that this is a 
slightly stronger rejection than the error bars quoted in the 
first paragraph of this Section, which correspond to a 1.7a 
rejection. The discrepancy is likely because renormalising 
the reduced statistic to 1 is not an entirely fair represen¬ 
tation of the error. If this renormalisation is not carried out, 
the rejection of the 600 km s“^ model rises to 1.9a. 

Wu et al. (1996) have pointed out that the inclusion 
of a uniform plane in modelling a dark matter halo may be 
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considered a good reflection of the lensing influence of large- 
scale structure. Following the model of Groom & Shanks 
(1999), we have considered the effects of including such a 
plane in our SIS model. Including a plane of dark matter 
with our SIS has no influence on the best-fit model or it’s 
error-bars. In fact, the only effect worth remarking is to 
slightly lower our rejection of low velocity dispersion models. 
The best-fit SIS model rejects 600 km s“^ at a 1.7a level. 
When a plane is included, this rejection drops to 1.5a. If 
we consider the group centre auto-correlation function (wcc) 
of Stevenson, Fong and Shanks (1988), then we might ex¬ 
pect the clustering of groups to have little effect on our sig¬ 
nal. Integrating under cUcc out to 10 arcmin suggests only a 
25 per cent contribution to our signal from the clustering of 
groups. 

There is still some debate over a number of the pa¬ 
rameters used in the modelling process. Changing the QSO 
number-count faint-end slope to /I = 0.34 raises the SIS 
model estimate of the velocity dispersion by ~ 10 per cent 
and lowering the index to = 0.24 lowers the estimate about 
10 per cent. Changing to a Concordance cosmology with 
Ilm = 0.3 and A = 0.7, to model lens-source separations 
raises the estimate by ~ 5 per cent. In the case of the NFW 
profile, lowering (5 to 0.24 lowers the mass estimate by ~ 
30 per cent, raising (5 to 0.34 raises the mass estimate by ~ 
50 per cent and changing to a Concordance cosmology raises 
the mass estimate by ~ 5 per cent. In short, cosmology does 
not really affect our estimates but the exact gradient of the 
faint-end slope of the QSO number-magnitude counts may 
well be important, especially for the NFW profile. Reducing 
the QSO number-count faint-end slope to = 0.15, a strong 
minimum value (being the differential counts slope deter¬ 
mined in Section 2.1) would reduce the NFW mass estimate 
by 70 per cent, bringing the NFW model mass prediction in 
line with a ACDM cosmology. It is possible that the magnifi¬ 
cation values generated by groups are large enough to reach 
so faintly into the QSO counts that the value of P is as low as 
the differential counts value. This could conceivably allow a 
ACDM cosmology combined with an NFW profile to explain 
the anti-correlation between QSOs and galaxy groups. 


5.3 Estimating Q.m in Groups 

In the Appendix, we also note how estimates of the average 
group mass, or velocity dispersion, can be turned into esti¬ 
mates of the mass density of the Universe that is associated 
with groups (see Croom & Shanks 1999). The sky density of 
our groups is around 0.85 deg“^ in both the NGC and SGC 
strips. Note that this is significantly larger than the density 
of Abell clusters (~ O.ldeg”^). Croom & Shanks (1999) 
have estimated the typical space density of such groups as 
3±1 X 10“^/i®Mpc“®. Using this value. Equation (A22) may 
be written 


Ilm — ( 


1125kms“ 


(7) 


assuming that the extent of the anti-correlation signal 
is 0 = 10 ± 2 arcmin (r = 1 ± Q.2h~^ Mpc). Equation (A23) 
is equivalently 



6 (arcmin) 

Figure 15. The cross-correlation between 2QZ QSOs and ob¬ 
jectively derived galaxy groups of limiting magnitude B = 20.5 
combined for both the southern and northern 2QZ strips together 
with best-fit models. Errors are standard deviations on the best¬ 
fitting NFW model derived from 500 mock QSO samples. 


^ _ AInfw ^ All.5 .c,., 

™ ~ 9.2 X IO 14 / 1 - 1 M 0 - 9.2 X IO 14 / 1 - 1 M 0 ^ ’ 

Using our best-fit estimates for the SIS (including scales 
between 10 and 40 arcsec) returns a value of Qm = 1.06±o;gi 
and for the NFW ~ 1.3, with large error. 

The large error in Croom & Shanks (1999) value for 
the space density of groups remains a dominant systematic 
in our estimates of and needs to be reassessed when 
group redshifts become available. Some additional error may 
be introduced by a lack of accurate redshift determinations 
for our galaxy groups. Groups of galaxies that are actu¬ 
ally greatly separated in redshift may accidentally align and 
thus be counted as a single halo, although it is unclear to 
what extent this contamination could influence our lensing 
results, as any associations along the line of sight still trace 
an increase in the intervening mass distribution. Certainly, 
associations of galaxies that are separated greatly in redshift 
will not have dark matter profiles like the SIS and NFW pro¬ 
files used in our modelling, being more like filaments than 
single haloes. Ray tracing of high resolution, N-body sim¬ 
ulations of the foreground mass distribution, where we can 
also apply our group detection algorithms, are needed to 
test the size of the anti-correlation expected under a specific 
cosmological model such as ACDM. 

Many of the inaccuracies in our method could be re¬ 
solved by a large sample of complementary galaxies and 
QSOs with redshift determinations and as such we await the 
publishing of the 2dF Galaxy Redshift Survey (Colless 2001) 
in its entirety. In the meantime, our estimates of Qm associ¬ 
ated with groups continue to appear as high as those found 
by Croom & Shanks. However, the errors have increased 
meaning the rejection of ~ 0.3 is only at the l-2a level 
and we have noted that there may be further systematic er- 
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Figure 16. The distribution of the strength of correlation for 
cross-correlations between 250 mock QSO catalogues and objec¬ 
tively derived galaxy groups. The mock QSO catalogues have been 
weighted to reflect a number of different possible SIS halo mod¬ 
els. The significance of the cross-correlation is measured out to 
10 arcmin for each mock QSO vs galaxy group cross-correlation 
function. Cross-correlating 2QZ QSOs against galaxy groups re¬ 
sulted in an anti-correlation of —2.9a', which is marked as the 
dashed line. 

rors still to be taken into account. Nevertheless, it is still 
tantalising that a a = 600 km s“^ dispersion is acceptable 
in at most 13/250 simulations for groups as numerous as 
those used here and this clearly motivates the application of 
this technique in larger QSO and galaxy group surveys. 

5.4 Discussion 

The approach to the modelling of QSO lensing adopted 
in this paper is different to that of some other authors. 
Williams & Irwin (1998) used the galaxy-QSO cross¬ 
correlation function coupled with the galaxy-galaxy auto¬ 
correlation function to derive constraints on Qm/ b, where b 
is the galaxy bias parameter and found their phenomeno¬ 
logical model largely agreed with sophisticated analytic 
approaches (Dolag & Bartelmann 1997; Sanz, Martinez- 
Gonzalez & Benitez 1997). In a similar study utilising 
the cluster-QSO cross-correlation function coupled with the 
cluster-cluster auto-correlation, Groom & Shanks (1999) ob¬ 
tained an estimate of Q/b ~ 3 —4. However, this estimate as¬ 
sumes that the total contribution to the oJcq anti-correlation 
comes from lensing by other clusters at a distance from the 
considered central cluster and that the contribution from the 
central cluster is negligible, which may not be true on scales 
a few arcminutes away from the cluster centre. The same 
criticism does not apply to the Williams and Irwin (1998) re¬ 
sult, it being detected on scales of tens of arcminutes, where 
galaxy-galaxy clustering will dominate uigq rather than the 
central galaxy halo. Their estimate of flm/b ~ 4-5 there¬ 
fore stands, although their statistical errors may still allow 
a value of Qm/6 ~ 1. 


Above, we determine the anti-correlation between 
galaxy groups and QSOs is best fit by a high-mass NFW 
model ( ~ 1.2 X 10^® h~^ Mq) or a high-velocity disper¬ 
sion SIS model (~ 1150 km s“^). Although the preference 
for the 1140 km s“^ SIS group velocity dispersion is only 
at the ~ 2a level in the current data, taking this result 
together with the previous independently derived result of 
BFS88 and also with the strong positive correlations seen 
in the brighter LBQS sample, it is clearly worth considering 
the implications if the amplitude of the anti-correlation was 
correct and caused wholly by weak lensing. Dynamical anal¬ 
ysis of 2dF galaxy z-space distortions, results in an estimate 
of /3 = /b ~ 0.43±0.07 (Peacock et al. 2001) and so mea¬ 

surement of Qm ~ 1 would imply b ~ 2.5. Although there is 
no immediate contradiction with the result presented a con¬ 
tradiction does arise when current CMB constraints on the 
mass power spectrum are included. These suggest that the 
galaxy power spectrum is approximately unbiased, implying 
Qm ~ 0.3, in contradiction at the 2a level with our best-fit 
result. 

There are other constraints on the mass of galaxy 
groups which are in contradiction with our best-fit velocity 
dispersion. In particular, Hoekstra et al. (2001) have used 
shear to measure weak lensing of galaxies behind GNOC 
groups, finding a ~ 300 km s“^ and M/L ^ 200hMQ/LQ. 
These translate into a value of Qm = 0.19. It is not clear 
why there is a difference between their results and ours. If 
the result of Hoekstra is correct, we would have to appeal 
to moderate statistical fluctuations to explain our high anti¬ 
correlation amplitude in an Qm = 0.3 model. Glearly shear 
studies behind nearby galaxy groups which also have QSO 
lensing data will be valuable. 

It should also be noted that the number of galaxies de¬ 
tected in N > 7 groups corresponds to only 7.1 per cent of 
the total number of galaxies. These exist, on average, within 
2.5 arcmin of the cluster centre. As mentioned in Section 3, 
our group centre against galaxy correlation function (wcg) 
agrees with Stevenson, Fong and Shanks (1988). Integrating 
ujcg to a radius of r < 10/ the extent of our anti-correlation 
signal, suggests that the total number of galaxies associated 
with our groups is actually ~ 20 per cent. If we assume that 
the M/L ratio of galaxies in clusters is the same as outside 
the groups then our estimate of Qm would rise by a factor of 
~ 5. However, in the ACDM model, for example, the galaxy 
distribution is expected to be anti-biased on small scales ie 
the M/L in clusters is expected to be higher by as much as 
~ 2 (Colin et al. 1999) and this introduces a further uncer¬ 
tainty into our best estimate of Qm. This suggests that our 
best estimate of Qm ~ 1 might be considered a weak lower 
limit and our l-2a value of Qm ~ 0.3 a strong lower limit. 
Any stronger conclusions from the present dataset await a 
more detailed test in a ACDM N-body simulation of the 
foreground mass distribution. 


6 CONCLUSIONS 

We have sought the effects of weak gravitational lensing in 
correlations between 2QZ QSOs and galaxy groups derived 
from the SDSS or from the Southern APM. We confirm 
that there is a distinct (~ 30-) anti-correlation between ob¬ 
jectively determined galaxy groups and QSOs (Boyle, Fong 
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& Shanks 1988). The anti-correlation is fit well by suppos¬ 
ing its cause is gravitational lensing through dark matter 
haloes, either NFW haloes or Singular Isothermal Spheres, 
but requires more mass than models with Q.m = 0.3 would 
suggest, at the ~ 1-2g level. Larger QSO samples could not 
only better distinguish the amplitudes of the anti-correlation 
predicted by cosmological models but may also be able to 
distinguish between different forms for dark matter halo pro¬ 
files. 

The observed anti-correlation does not appear to be 
caused by a selection effect due to the limited spacing of 2dF 
fibres. We also rule out the idea that intervening dust causes 
the signal (Boyle, Fong & Shanks 1988). Our 95 per cent 
upper limits on reddening in bj — r are 0.012, which cor¬ 
responds to Ab < 0.04 mag assuming the usual Galactic 
dust law. To explain the anti-correlation by dust would need 
Ab ~ 0.2 mag (see BFS88); the dust hypothesis could then 
only be saved by appealing to either a smooth distribution of 
grey dust or a patchy distribution of heavily line of sight dis¬ 
tributed dust around galaxy groups. It is also not straight¬ 
forward to rule out the hypothesis that both dust and lensing 
play some combined role in producing the anti-correlation 
signal. 

It seems that weak lensing remains the likely expla¬ 
nation for the anti-correlation between QSOs and galaxy 
groups. The strength of the anti-correlation suggests that 
there may be more mass hidden in galaxy groups than many 
previous estimates require, making further study of this phe¬ 
nomenon worthwhile. The models used in Section 5.2 to de¬ 
scribe the form of the correlation function remain simplis¬ 
tic averages and it would be worth performing some large, 
high-resolution simulations of the expected lensing influence 
of foreground mass on QSO distributions in different cos¬ 
mologies. Accurate ray tracing through N-body simulations 
would predict the expected anti-correlation for different cos¬ 
mologies. Additionally, it is probable that a large proportion 
of our groups are really chance alignments of galaxies that 
are actually greatly separated in redshift. Running our group 
detection procedure using galaxies from a large simulation, 
should allow us to determine how many of our groups are 
genuine haloes and how many are associations or filamen¬ 
tary in structure. This, in turn, would allow us to produce 
a more accurate model of the mass profile we are fitting to 
the data. 

On the observational side, new faint QSO surveys linked 
to the 2QZ should definitively determine the faint-end slope 
of the deep QSO number-magnitude relation. Redshift in¬ 
formation in the 2dF Galaxy Redshift Survey should im¬ 
prove our estimate of flm a great deal. Using forthcoming 
catalogues of 2dFGRS galaxy groups (Eke et ah, in prepa¬ 
ration) it should be possible to accurately look at the anti¬ 
correlation amplitude as a function of group/cluster rich¬ 
ness to try and further distinguish the masses associated 
with groups and clusters. Groups determined from galax¬ 
ies of known redshift have the added advantage of tracing 
definite dark matter haloes, rather than filaments or allign- 
ments. Certainly, more statistical power is also needed, mo¬ 
tivating extending the 2dF QSO over a wider area and to 
greater depth. The large area of the Sloan Digital Sky Sur¬ 
vey, containing both QSO and galaxy samples, may be use¬ 
ful for both accurate grouping of galaxies in redshift space 
and measuring the extent of the anti-correlation on small 


scales, where the mass and its form is best constrained. No¬ 
tably, the SDSS will contain significant samples of QSOs 
brighter than bj ~ 19, allowing us to look for confirmation 
of the expected positive correlation between galaxies and 
QSOs brighter than the knee of the QSO number counts. 
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APPENDIX A: STATISTICAL LENSING OF 
QSOS THROUGH DARK MATTER HALOES 

A1 General lensing 

As Einstein noted (Einstein 1915), a mass M will bend a 
ray of light passing at impact parameter b through an angle 
a (Schneider, Ehlers & Falco 1992) 


_ 4GM(< b) _ D, 
bc^ ^ As 


(Al) 


where Di,Ds and Dis are the angular diameter distances of 
the lens as measured by the observer, the source as measured 
by the observer and the source as measured by the lens, 
respectively; 6 is the angle between the image, the observer 
and the centre of the lens; 6q is the angle between the source, 
the observer and the centre of the lens. Flux conservation 
implies each image will be magnified by a factor 


n = 


e de 


(A2) 


(Turner, Ostriker & Gott 1984). Lensing influences a sam¬ 
ple of background objects in two competing ways. Fainter 
objects are lensed into a magnitude limitied sample, increas¬ 
ing the number density of that sample but the area behind 
the lens is proportionately expanded, reducing the sample’s 


number density. Narayan (1989) quantified this effect as a 
‘net enhancement factor’ 


1 N{< m J- 2.5log{fj,)) 
/r N{< m) 


(A3) 


The Number-magnitude relation can be approximated as a 
power law, with N{< m) oc 10^"“ (Boyle, Fong & Shanks 
1988), allowing us to express the enhancement factor as 


<7 


1 


j^Q/3(m+2.5iog((i)) 


= M 


2.5/3-1 


(A4) 


(Groom 1997). Now, the net enhancement factor is the ratio 
of the observed (lensed) flux and the true (unlensed) flux. 
The correlation function ui{9) may be expressed as the ratio 
of observed pairs of objects to expected pairs of objects. 
Typically, 1 is subtracted from the correlation function to 
account for the expected normal background of pairs. Hence: 


io{9) = q — 1 — ^ — 1 (A-5) 

Equation (A5) dictates u!{9) = 0 when /? = 0.4. For 
higher values of (5 we would observe a correlation, and for 
lower values, an anti-correlation. Thus, the lensing effect is 
dependent on the slope of the number-magnitude relation. 


A2 Dark matter profiles 

It is simple to integrate a Singular Isothermal Sphere (SIS) 
prohle out to an impact parameter b (Groom 1997) to de¬ 
termine that it will magnify background sources by a factor 


n = 


9 




(A6) 


with c being the speed of light and a the velocity disper¬ 
sion of the SIS. The term 47r^(-)^ in the denominator is 
the Einstein Radius. The Einstein Radius essentially marks 
the boundary between the weak and strong lensing regimes. 
Equations (A5) and (A6) can be combined to constrain the 
predicted form of the correlation function from magnifica¬ 
tion through an SIS prohle. 

N-body simulations have provided a universal (NEW) 
density prohle for dark matter haloes (Navarro, Frenk & 
White 1995, 1996, 1997) that has been independently obser- 
vationally conhrmed with some success (Bartelmann et al. 
1998; Thomas et al. 1998). We have also constructed a sim¬ 
ple model to determine the form of the correlation function 
based on lensing through such a halo. 

The NEW prohle in the form 


p(r) = 


be pc 


21(1 + 21 ) 

rs5 ^ rs5 z 


(A7) 


where rg is a representative radial scale and pc is the 
critical density or the universe at the redshift of the dark 
matter, appears to be a reasonable description of haloes 
spanning 9 orders of magnitude in mass, from globular clus¬ 
ters to galaxy clusters. Ideally, we would wish to study sta¬ 
tistical lensing utilising just such a realistic density prohle. 
Following Bartelmann (1996), we rewrite the prohle as 
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p{x) = 


Ps 

*(1 + xY 


(A8) 


and consider lensing around the profile at a radius of 
the impact parameter b, meaning that x is defined by 


A3 Cosmology: Determining fim in Groups 

We will now briefly outline how the mass density of the 
Universe might be determined from our study, noting that 
the value thus derived is purely a measure of flm in galaxy 
groups. Q,m may be defined 


b 


* Ta 


Die 

Xs 


(A9) 


where D\ and 9 are defined as for the SIS, above. One 
can then combine Equations (A2) and (A9) to yield: 


p. = 


Vs .2 X d® 

A ^d^ 


(AlO) 


Now, Maoz et al. (1997) determine the characteristic 
scale Xa from an empirical fit to Figure 9 of Navarro et al. 
(1997) as 


where h is the Hubble constant in units of 
100 km s“^Mpc“^ and 7 varies between cosmogonies (for 
CDM, 7 ~ i). Maoz et al. (1997) go on to rewrite Equation 
(Al) as 




po _ SttG 


(A21) 


where po can be estimated as the product of the space den¬ 
sity of galaxy groups and the average mass of a dark matter 
halo. The sky density of galaxy groups can be measured. 
Equations (A5) and (A 6 ) relate the average velocity disper¬ 
sion of the SIS to the observed correlation function. The 
average mass of the SIS may be derived from the average 
velocity dispersion, yielding an estimate of Qm in groups of 


„ SttG , . Stt^ 2 / a 

(A22) 

where n is the space density of galaxy groups and r refers to 
the extent of the dark matter halo or, roughly, the extent of 
correlation. Similarly, we can express the correlation func¬ 
tion in terms of the mass within 1.5h“^ Mpc of the centre of 
an NFW halo. Theoretically, then, Qrn is alternatively given 
by 


4 GMi .5 g{x) 

C^XagA-^Mpc/ra) X 


(A12) 

(A13) 


where Mi .5 is the mass within 1.5/i ^ Mpc of the centre 
of the halo and g{x) is given by Bartelmann (1996) as: 


„ SttG , SttG 

Urn = ^^^nMNFW > 


^0 


(A23) 


where Mi ,5 is the mass within 1.5h ^ Mpc of the centre 
of an NFW halo. 


/ \ X 2 _i 

g{x) = In - + ytan 

/ \ n ® 2 , _1 

g{x) = In - + 


X — 1 
a; -I- 1 

1 — X 
1 -I- a: 


(x > 1 ) 


(x < 1 ) 


g{x) = In - -f 1 (a; = 1 ) 


(A14) 

(A15) 

(A16) 


Gombining Equations (Al) and (A12), we can deduce 
As 4 GMi,5 g{x) 


Q — 0 _ 

A rsC^p(1.5Mpc/rs) x 


dSq d^ Dla 


4GMi,5 


da; da; A raAg{1.5Mpc/ra) da; 

and substitute in Equation (A9) to ultimately find 


(A17) 

(A18) 


^ ^ _ As 4 GMi.5 g{x) 
X A A rac'^g{1.5Mpc/ra) x^ 

^ = II _ As 4 GMi.5 

da; A A raC^gil.bMpc/rs) da; 


(A19) 

(A20) 


which can be readily substituted into Equation (AlO) 
and numerically solved to derive the magnification. The 
magnification dictates the expected form of the correlation 
function as outlined in Equation (A5). 
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